home *** CD-ROM | disk | FTP | other *** search
/ Sky at Night 2006 August / SAN CD 8-2006 CD-ROM 15.iso / pc / Media / Bonus Content / Jupiter's Moons / jupiter.js < prev    next >
Encoding:
Text File  |  2006-01-26  |  14.5 KB  |  494 lines
  1. /*
  2. Computational engine, setup and DHTML graphics:
  3. copyright (c) Adrian R. Ashford, 9th January 2006
  4.  
  5. JavaScript(s) used with permission by BBC Sky at Night magazine.
  6.  
  7. DHTML date & time validation scripts:
  8. Sandeep Tamhankar
  9. stamhankar@hotmail.com
  10. http://javascript.internet.com
  11. */
  12.  
  13. x_origin = 318;
  14. y_origin = 134;
  15. scale_factor = 12;
  16. o_x = 1;
  17. o_y = 1;
  18.  
  19. function setup()
  20. {
  21. var nowdate = new Date();
  22. var utc_day = nowdate.getUTCDate();
  23. var utc_month = nowdate.getUTCMonth() + 1;
  24. var utc_year = nowdate.getUTCFullYear();
  25. zone = nowdate.getTimezoneOffset() / 1440;
  26. var utc_hours = nowdate.getUTCHours();
  27. var utc_mins = nowdate.getUTCMinutes();
  28. var utc_secs = nowdate.getUTCSeconds();
  29. utc_mins += utc_secs / 60.0;
  30. utc_mins = Math.floor((utc_mins + 0.5));
  31.  
  32. if (utc_mins > 59) utc_mins = 59;
  33.  
  34. //document.redspot.date_txt.value = utc_day + "/" + utc_month + "/" + utc_year;
  35. //document.redspot.ut_h_m.value = utc_hours + ":" + utc_mins
  36.  
  37. document.redspot.m1.value = doubleDigit(utc_day);
  38. document.redspot.m2.value = doubleDigit(utc_month);
  39. document.redspot.m3.value = utc_year;
  40. document.redspot.m4.value = doubleDigit(utc_hours);
  41. document.redspot.m5.value = doubleDigit(utc_mins);
  42.  
  43. jupiter(0,0);
  44. }
  45.  
  46. function doubleDigit(iNum){
  47.     if(iNum < 10)
  48.         return '0' + iNum;    
  49.     return iNum;    
  50. }
  51.  
  52. function proper_ang(big)
  53. {
  54. with (Math)
  55. {
  56. var tmp = 0;
  57. if (big > 0)
  58. {
  59. tmp = big / 360.0;
  60. tmp = (tmp - floor(tmp)) * 360.0;
  61. }
  62. else
  63. {
  64. tmp = ceil(abs(big / 360.0));
  65. tmp = big + tmp * 360.0;
  66. }
  67. }
  68. return tmp;
  69. }
  70.  
  71. function round_10(num)
  72. {
  73. return Math.floor((num + 0.05) * 10) / 10;
  74. }
  75.  
  76. function isValidDate()
  77. {
  78.     dateStr = document.redspot.m1.value+ "/" + document.redspot.m2.value +"/"+ document.redspot.m3.value;
  79.     timeStr = document.redspot.m4.value +":"+ document.redspot.m5.value;
  80.     
  81. if (IsValidTime(timeStr) == true)
  82. {
  83. var datePat = /^(\d{1,2})(\/|-)(\d{1,2})\2(\d{4})$/;
  84. var matchArray = dateStr.match(datePat);
  85. if (matchArray == null) { 
  86. alert("Date is not in a valid format.") 
  87. return false; 
  88. } 
  89. var month = matchArray[3];
  90. var day = matchArray[1];
  91. var year = matchArray[4];
  92. if (month < 1 || month > 12) {
  93. alert("Month must be between 1 and 12.");
  94. return false; 
  95. } 
  96. if (year < 2005 || year > 2015) {
  97. alert("Year must be between 2005 and 2015.");
  98. return false; 
  99. } 
  100. if (day < 1 || day > 31) {
  101. alert("Day must be between 1 and 31."); 
  102. return false; 
  103. } 
  104. if ((month == 4 || month == 6 || month == 9 || month == 11) && day == 31) {
  105. alert("Month " + month + " doesn't have 31 days!") 
  106. return false 
  107. } 
  108. if (month == 2) {
  109. var isleap = (year % 4 == 0 && (year % 100 != 0 || year % 400 == 0)); 
  110. if (day > 29 || (day == 29 && !isleap))
  111. { 
  112. alert("February " + year + " doesn't have " + day + " days!"); 
  113. return false; 
  114. } 
  115. } 
  116. if (month < 10 && month.length == 1) month = "0" + month;
  117. if (day < 10 && day.length == 1) day = "0" + day;
  118.  
  119. document.redspot.m1.value = day;
  120. document.redspot.m2.value = month;
  121. document.redspot.m3.value = year;
  122. var dt_str = day +"/"+ month + "/" + year;
  123. if ((dt_str.substring(2,3) != "/") || (dt_str.substring(5,6) != "/"))
  124. {
  125. alert ("Date is not in a valid format.");
  126. return false;
  127. }
  128. jupiter(0,0);
  129. return true;
  130. }
  131. else
  132. {
  133. return false;
  134. }
  135. }
  136.  
  137. function IsValidTime(timeStr)
  138. {
  139. var timePat = /^(\d{1,2}):(\d{2})(:(\d{2}))?(\s?(AM|am|PM|pm))?$/;
  140. var matchArray = timeStr.match(timePat);
  141. if (matchArray == null)
  142. {
  143. alert("Time is not in a valid format.");
  144. return false;
  145. }
  146. var hour = matchArray[1];
  147. var minute = matchArray[2];
  148. if (hour < 0 || hour > 23)
  149. {
  150. alert("Hour must be between 0 and 23.");
  151. return false;
  152. }
  153. if (minute < 0 || minute > 59)
  154. {
  155. alert ("Minute must be between 0 and 59.");
  156. return false;
  157. }
  158. document.redspot.m4.value = hour;
  159. document.redspot.m5.value = minute;
  160. var tm_str =hour + ":" + minute;
  161. if ((tm_str.substring(2,3) != ":") && (dt_str.length != 5))
  162. {
  163. alert ("Time is not in a valid format.");
  164. return false;
  165. }
  166. return true;
  167. }
  168.  
  169. function isBinocular()
  170. {
  171. o_x = 1;
  172. o_y = 1;
  173. jupiter(0,0);
  174. return true;
  175. }
  176.  
  177. function isReflector()
  178. {
  179. o_x = -1;
  180. o_y = -1;
  181. jupiter(0,0);
  182. return true;
  183. }
  184.  
  185. function isRefractor()
  186. {
  187. o_x = -1;
  188. o_y = 1;
  189. jupiter(0,0);
  190. return true;
  191. }
  192.  
  193. /*
  194. function time_change(tmp)
  195. {
  196. if (isValidDate(document.redspot.date_txt.value) == true)
  197. {
  198. var jd_temp, zz, ff, alpha, aa, bb, cc, dd, ee;
  199. var calendar_day, calendar_month, calendar_year;
  200. var int_day, hours, minutes;
  201.  
  202. var tm_as_str, ut_hrs, ut_mns, frac_day;
  203. var jd = julian_date();
  204.  
  205. ut_hrs = document.redspot.m4.value
  206. ut_mns = document.redspot.m5.value
  207. frac_day = ut_hrs / 24.0 + ut_mns / 1440.0;
  208.  
  209. with (Math) {
  210.  
  211. jd_temp = jd + frac_day + tmp / 24.0 + 0.5;
  212.  
  213. zz = floor(jd_temp);
  214. ff = jd_temp - zz;
  215. alpha = floor((zz - 1867216.25) / 36524.25);
  216. aa = zz + 1 + alpha - floor(alpha / 4);
  217. bb = aa + 1524;
  218. cc = floor((bb - 122.1) / 365.25);
  219. dd = floor(365.25 * cc);
  220. ee = floor((bb - dd) / 30.6001);
  221. calendar_day = bb - dd - floor(30.6001 * ee) + ff;
  222. calendar_month = ee;
  223. if (ee < 13.5) calendar_month = ee - 1;
  224. if (ee > 13.5) calendar_month = ee - 13;
  225. calendar_year = cc;
  226. if (calendar_month > 2.5) calendar_year = cc - 4716;
  227. if (calendar_month < 2.5) calendar_year = cc - 4715;
  228. int_day = floor(calendar_day);
  229. hours = (calendar_day - int_day) * 24;
  230. minutes = floor((hours - floor(hours)) * 60 + 0.5);
  231. hours = floor(hours);
  232. if (minutes > 59)
  233. {minutes = 0; hours = hours + 1;}
  234. if (calendar_month < 10) calendar_month = "0" + calendar_month;
  235. if (int_day < 10) int_day = "0" + int_day;
  236. if (hours < 10) hours = "0" + floor(hours);
  237. if (minutes < 10) minutes = "0" + minutes;
  238. }
  239. jupiter();
  240. return true;
  241. }
  242. else
  243. {
  244. return false;
  245. }
  246. }
  247. */
  248. function julian_date()
  249. {
  250. var dt_as_str, mm, dd, yy;
  251. var yyy, mmm, a, b;
  252.  
  253. mm = Number(document.redspot.m2.value);
  254. dd = Number(document.redspot.m1.value);
  255. yy = Number(document.redspot.m3.value);
  256.  
  257. with (Math) {
  258. var yyy=yy;
  259. var mmm=mm;
  260. if (mm < 3)
  261. {
  262. yyy = yy - 1;
  263. mmm = mm + 12;
  264. }
  265. a = floor(yyy/100);
  266. b = 2 - a + floor(a/4);
  267.  
  268. return floor(365.25*yyy) + floor(30.6001*(mmm+1)) + dd + 1720994.5 + b;
  269. }
  270. }
  271.  
  272. function jupiter(oMenID, iVal)
  273. {
  274. /*document.redspot.eval(oMenID).value = doubleDigit(Number(document.redspot.eval(oMenID).value) + iVal)*/
  275. switch(oMenID){
  276.     case 'm1':
  277.         document.redspot.m1.value = doubleDigit(Number(document.redspot.m1.value) + iVal);
  278.         break;
  279.     case 'm2':
  280.         document.redspot.m2.value = doubleDigit(Number(document.redspot.m2.value) + iVal);
  281.         break;    
  282.     case 'm3':
  283.         document.redspot.m3.value = doubleDigit(Number(document.redspot.m3.value) + iVal);
  284.         break;    
  285.     case 'm4':
  286.         document.redspot.m4.value = doubleDigit(Number(document.redspot.m4.value) + iVal);
  287.         break;    
  288.     case 'm5':
  289.         document.redspot.m5.value = doubleDigit(Number(document.redspot.m5.value) + iVal);
  290.         break;
  291. }
  292.  
  293. var i,trans,io_layer,europa_layer,ganymede_layer,callisto_layer;
  294. var tm_as_str, ut_hrs, ut_mns, frac_day;
  295. var moon_separation, shadow_separation;
  296. var factor = 1.071374;
  297.  
  298. var jd = julian_date();
  299.  
  300. ut_hrs = document.redspot.m4.value
  301. ut_mns = document.redspot.m5.value
  302. frac_day = ut_hrs / 24.0 + ut_mns / 1440.0;
  303.  
  304. with (Math) {
  305.  
  306. var RAD = 180.0/PI;
  307.  
  308. var day_start = jd + 1.1 / 1440 - 2451545.0;
  309. var days = day_start + frac_day;
  310.  
  311. var V = proper_ang(172.74 + 0.00111588 * days);
  312. V = 0.329 * sin(V / RAD);
  313. var M = proper_ang(357.529 + 0.9856003 * days) / RAD;
  314. var N = proper_ang(20.02 + 0.0830853 * days + V) / RAD;
  315. var J = proper_ang(66.115 + 0.9025179 * days - V);
  316. var A = 1.915 * sin(M) + 0.02 * sin(2 * M);
  317. var B = 5.555 * sin(N) + 0.168 * sin(2 * N);
  318. var K = (J + A - B) / RAD;
  319. var R = 1.00014 - 0.01671 * cos(M) - 0.00014 * cos(2 * M);
  320. var r = 5.20872 - 0.25208 * cos(N) - 0.00611 * cos(2 * N);
  321. var e_to_j = sqrt(pow(r,2) + pow(R,2) - 2 * r * R * cos(K));
  322. document.redspot.jupiter_dist.value = floor((e_to_j + 0.005) * 100.0) / 100.0;
  323.  
  324. var ang_jove = 0.0009543 / e_to_j * RAD * 30.0;
  325. document.redspot.jupiter_size.value = round_10(ang_jove * 120.0);
  326. var pha = asin(R / e_to_j * sin(K));
  327. var lambda = proper_ang(34.35 + 0.083091 * days + V + B);
  328. var D_s = 3.12 * sin ((lambda + 42.8) / RAD);
  329. var D_e = D_s - 2.22 * sin(pha) * cos((lambda + 22) / RAD) - 1.3 * (r - e_to_j) / e_to_j * sin((lambda - 100.5) / RAD);
  330. var diff = sin((D_e - D_s) / RAD);
  331. D_e = sin(D_e / RAD);
  332. D_s = sin(D_s / RAD);
  333. pha = pha * RAD;
  334. var m_txt = new Array (4);
  335. var moons_x = new Array (4);
  336. var moons_y = new Array (4);
  337. var shadow_x = new Array (4);
  338. var shadow_y = new Array (4);
  339. var moon_shadow_x = new Array (4);
  340. var moon_shadow_y = new Array (4);
  341. var u = new Array (4);
  342. var u1_u = proper_ang(163.8069 + 203.4058646 * (days - e_to_j/173) + pha - B);
  343. var u2_u = proper_ang(358.414 + 101.2916335 * (days - e_to_j/173) + pha - B);
  344. var u3_u = proper_ang(5.7176 + 50.234518 * (days - e_to_j/173) + pha - B);
  345. var u4_u = proper_ang(224.8092 + 21.48798 * (days - e_to_j/173) + pha - B);
  346. var G = proper_ang(331.18 + 50.310482 * (days - e_to_j/173)) / RAD;
  347. var H = proper_ang(87.45 + 21.569231 * (days - e_to_j/173)) / RAD;
  348. var sh_factor_x = sin(pha / RAD);
  349. var sh_factor_y = cos(pha / RAD);
  350. u[0] = (u1_u + 0.473 * sin(2 * (u1_u - u2_u) / RAD)) / RAD;
  351. u[1] = (u2_u + 1.065 * sin(2 * (u2_u - u3_u) / RAD)) / RAD;
  352. u[2] = (u3_u + 0.165 * sin(G)) / RAD;
  353. u[3] = (u4_u + 0.843 * sin(H)) / RAD;
  354. var r1 = 5.9057 - 0.0244 * cos(2 * (u1_u - u2_u) / RAD);
  355. var r2 = 9.3966 - 0.0882 * cos(2 * (u2_u - u3_u) / RAD);
  356. var r3 = 14.9883 - 0.0216 * cos(G);
  357. var r4 = 26.3627 - 0.1939 * cos(H);
  358.  
  359. moons_x[0] = r1 * sin(u[0]);
  360. moons_x[1] = r2 * sin(u[1]);
  361. moons_x[2] = r3 * sin(u[2]);
  362. moons_x[3] = r4 * sin(u[3]);
  363. moons_y[0] = -(r1 * cos(u[0]) * D_e) * factor;
  364. moons_y[1] = -(r2 * cos(u[1]) * D_e) * factor;
  365. moons_y[2] = -(r3 * cos(u[2]) * D_e) * factor;
  366. moons_y[3] = -(r4 * cos(u[3]) * D_e) * factor;
  367. shadow_x[0] = -(r1 * sh_factor_x);
  368. shadow_x[1] = -(r2 * sh_factor_x);
  369. shadow_x[2] = -(r3 * sh_factor_x);
  370. shadow_x[3] = -(r4 * sh_factor_x);
  371. shadow_y[0] = r1 * sh_factor_y * diff * factor;
  372. shadow_y[1] = r2 * sh_factor_y * diff * factor;
  373. shadow_y[2] = r3 * sh_factor_y * diff * factor;
  374. shadow_y[3] = r4 * sh_factor_y * diff * factor;
  375. moon_shadow_x[0] = moons_x[0] + shadow_x[0];
  376. moon_shadow_y[0] = moons_y[0] + shadow_y[0];
  377. moon_shadow_x[1] = moons_x[1] + shadow_x[1];
  378. moon_shadow_y[1] = moons_y[1] + shadow_y[1];
  379. moon_shadow_x[2] = moons_x[2] + shadow_x[2];
  380. moon_shadow_y[2] = moons_y[2] + shadow_y[2];
  381. moon_shadow_x[3] = moons_x[3] + shadow_x[3];
  382. moon_shadow_y[3] = moons_y[3] + shadow_y[3];
  383.  
  384. var io_x = x_origin + moons_x[0] * scale_factor * o_x;
  385. var io_y = y_origin - moons_y[0] * scale_factor * o_y;
  386. var europa_x = x_origin + moons_x[1] * scale_factor * o_x;
  387. var europa_y = y_origin - moons_y[1] * scale_factor * o_y;
  388. var ganymede_x = x_origin + moons_x[2] * scale_factor * o_x;
  389. var ganymede_y = y_origin - moons_y[2] * scale_factor * o_y;
  390. var callisto_x = x_origin + moons_x[3] * scale_factor * o_x;
  391. var callisto_y = y_origin - moons_y[3] * scale_factor * o_y;
  392.  
  393. var io_shadow_x = x_origin + moon_shadow_x[0] * scale_factor * o_x;
  394. var io_shadow_y = y_origin - moon_shadow_y[0] * scale_factor * o_y;
  395. var europa_shadow_x = x_origin + moon_shadow_x[1] * scale_factor * o_x;
  396. var europa_shadow_y = y_origin - moon_shadow_y[1] * scale_factor * o_y;
  397. var ganymede_shadow_x = x_origin + moon_shadow_x[2] * scale_factor * o_x;
  398. var ganymede_shadow_y = y_origin - moon_shadow_y[2] * scale_factor * o_y;
  399. var callisto_shadow_x = x_origin + moon_shadow_x[3] * scale_factor * o_x;
  400. var callisto_shadow_y = y_origin - moon_shadow_y[3] * scale_factor * o_y;
  401.  
  402. io_layer = 12;
  403. europa_layer = 13;
  404. ganymede_layer = 14;
  405. callisto_layer = 15;
  406.  
  407. if (1.571 < u[0] && u[0] < 4.712) io_layer = 9;
  408. if (1.571 < u[1] && u[1] < 4.712) europa_layer = 7;
  409. if (1.571 < u[2] && u[2] < 4.712) ganymede_layer = 5;
  410. if (1.571 < u[3] && u[3] < 4.712) callisto_layer = 3;
  411.  
  412. for (i=0; i<4; i++)
  413. {
  414. m_txt[i] = " lies " + abs(round_10(moons_x[i] * ang_jove)) + " arcmin ";
  415. if (moons_x[i] > 0)
  416. {m_txt[i] = m_txt[i] + "west";}
  417. else
  418. {m_txt[i] = m_txt[i] + "east";}
  419. m_txt[i] = m_txt[i] + " and " + abs(floor(moons_y[i] * ang_jove * 60.0 + 0.5)) + " arcsec ";
  420. if (moons_y[i] > 0)
  421. {m_txt[i] = m_txt[i] + "north";}
  422. else
  423. {m_txt[i] = m_txt[i] + "south";}
  424. m_txt[i] = m_txt[i] + " of Jupiter's centre.";
  425.  
  426. shadow_separation = sqrt(pow((moons_x[i] - shadow_x[i]),2) + pow((moons_y[i] - shadow_y[i]),2));
  427. moon_separation = sqrt(pow(moons_x[i],2) + pow(moons_y[i],2));
  428. moon_shadow_separation = sqrt(pow(moon_shadow_x[i],2) + pow(moon_shadow_y[i],2));
  429.  
  430. if (shadow_separation < 1 && (u[i] > 1.58 && u[i] < 4.71))
  431. {
  432. m_txt[i] = " is currently in Jupiter's shadow.";
  433. if (i == 0) io_layer = 1;
  434. if (i == 1) europa_layer = 1;
  435. if (i == 2) ganymede_layer = 1;
  436. if (i == 3) callisto_layer = 1;
  437. }
  438.  
  439. if (moon_separation < 1 && (u[i] > 2.79 && u[i] < 3.49))
  440. m_txt[i] = " is currently occulted by Jupiter.";
  441. if (moon_separation < 1 && (u[i] > 5.93 || u[i] < 0.35))
  442. m_txt[i] = " is currently transiting Jupiter.";
  443.  
  444. if (moon_shadow_separation < 1 && (u[i] > 5.23 || u[i] < 1.04))
  445. m_txt[i] = " : SHADOW TRANSIT IN PROGRESS!";
  446. }
  447.  
  448. document.getElementById('io').style.left = io_x;
  449. document.getElementById('io').style.top = io_y;
  450. document.getElementById('io').style.zIndex = io_layer;
  451.  
  452. document.getElementById('europa').style.left = europa_x;
  453. document.getElementById('europa').style.top = europa_y;
  454. document.getElementById('europa').style.zIndex = europa_layer;
  455.  
  456. document.getElementById('ganymede').style.left = ganymede_x;
  457. document.getElementById('ganymede').style.top = ganymede_y;
  458. document.getElementById('ganymede').style.zIndex = ganymede_layer;
  459.  
  460. document.getElementById('callisto').style.left = callisto_x;
  461. document.getElementById('callisto').style.top = callisto_y;
  462. document.getElementById('callisto').style.zIndex = callisto_layer;
  463.  
  464. document.redspot.x1.value = m_txt[0];
  465. document.redspot.x2.value = m_txt[1];
  466. document.redspot.x3.value = m_txt[2];
  467. document.redspot.x4.value = m_txt[3];
  468.  
  469. var System_I = proper_ang(210.98 + 877.8169088 * (day_start - e_to_j/173) + pha - B);
  470. var System_II = proper_ang(187.23 + 870.1869088 * (day_start - e_to_j/173) + pha - B);
  471.  
  472. var Red_Spot = 104.0; // Red Spot's current longitude in degrees
  473.  
  474. var t_h = 0;
  475. var t_m = 0;
  476. var diff = proper_ang(Red_Spot - System_II);
  477. var rs_time = new Array(3);
  478. for (i=0; i<3; i++)
  479. {
  480. trans = diff / 870.1869 * 24.0;
  481. t_h = floor(trans);
  482. t_m = floor(60.0 * (trans - t_h) + 0.5);
  483. if (t_m > 59) t_m = 59;
  484. if (t_m < 10) t_m = "0" + t_m;
  485. if (t_h < 10) t_h = "0" + t_h;
  486. rs_time[i] = t_h + ":" + t_m + " UT";
  487. diff += 360.0;
  488. }
  489. document.redspot.rs_transit1.value = rs_time[0];
  490. document.redspot.rs_transit2.value = rs_time[1];
  491. document.redspot.rs_transit3.value = rs_time[2];
  492. if (t_h > 23) document.redspot.rs_transit3.value = ".........";
  493. }
  494. }